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INTRODUCTION. 



This publication was written for those persons who are interested in the 
application of wind energy and who want to know how to design the blade 
shape of a windmill rotor. We have received many requests on this issue 
from parties inside and outside the Netherlands. In writing this booklet 
we were uncertain as to the level at which it should be written. 
So some might find it too easy, others too difficult. 

We would like to emphasize that the design procedure, as 
given in one of the chapters, is very simple. It is important however that 
a number of basic ideas and concepts are well understood, before an attempt 
is made to design a rotor. Therefore quite a lot of attention is given to 
explaining lift, drag, rotor characteristics etc. So we ask the reader to 
be tolerant and patient. Those who are not familiar with the basic concepts 
we urge to go on and not to be afraid of the first formulas and graphs pre- 
sented here. 

Although a good rotor can be designed with the procedure as 
presented here, a few things should be noted. In the selection of a rotor 
type, in terms of design speed and radius, the load characteristics and 
wind availability must be taken into account. We have seen various good ro- 
tors coupled to wrong types of loads or to too high loads. Although it is 
possible to change the number of revolutions per minute with a transmission, 
this will not solve the problems that may arise when the selected design fi- 
gures are basically incorrect. A second remark is that rotors can be very 
dangerous. No strength calculations are given in this book, but remember that 
the centrifugal forces can make a rotor explode if it is not strong enough. 
Touching a rotor during operation will lead to serious injury. 
The availability of certain materials and technologies can be taken into ac- 
count in the earliest stages of the design. We therefore hope chat, with this 
book, the reader will be able to design a rotor that can be manufactured with 
the means and technologies as are locally available. 



LIST OF SYMBOLS 
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a constant 

n 2 

A area m 

B number of blades 

c chord m 

C, drag coefficient 
d 

lift coefficient 

Cp power coefficient 

C torque coefficient 
Q 

d diameter m 

D drag N 

E energy J 

E energy per volume Jm 
v 

f plate bending of arched steel plate m 

L lift N 

m mass kg 

n number of revolutions per second s ^ 

P power W 

Q torque Nm 

R rotor radius m 

r local radius m 

Re Reynolds number 

Re Re for B=C, =r=V=l 

n 1 03 _l 

u tangential blade speed at radius r ms 

-1 

V velocity, speed ms 

-1 

V undisturbed wmdspeed ms 

W relative velocity to rotor blade ms~* 

a angle of attack. 

a design value for angle of attack 

> blade angle, blade setting 

n factor for blade number effect on C 

B " 

A tip-speed ratio 

X design value for tip-speed ratio 
o 

A local speed ratio at radius r 

r .... 2-1 

v kinematic viscosity m s 

-3 

p density kgm 

3 - 1 

$> volume flow m s 

v 

$ angle between plane of rotation and relative 

flow velocity at the rotor blades 

ft rotor angular velocity s 
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UNITS OF MEASURE. 



The units used in this publication belong to, or are based on the so-cailed 
Systeme International d'Unites (SI). 

Those who are not familiar with these units or who have data given in units of 
other systems (for example windspeed in mph), we here give a short list with the 
conversion factors for the units that are most relevant for the design of windmill 
rotors . 



length 
area 
volume 
speed 

mass 

force 

torque 

energy 

power 



1 m 



i m 



* 3.28 ft 

= 10.76 ft' 



I m 3 =35.31 ft 3 »264.2 gallons 
I nis ' = 2.237 mph 
1 knot = 0.5144 ms = 1.15 mph 
1 kg = 2.205 lb 

= 0.225 lbf = 0. 102 kgf 
» 0.738 ft lbf 



1 N 
1 Nm 
1 J 
I W 
1 hp 
1 pk 



= 0.239 calories = 0.2777* 10 
-1 



-6 



kWh=lNm 



= 1 watt = 1 Js 
= 0.7457 kW 
= 0.7355 kW 



= 0.738 ft lbf s ^lNms 1 
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1. WIND ENERGY - WIND POWER. 

Wind is air in motion. The air has mass, but its density is low. 
When mass is moving with velocity V it has kinetic energy expressed 
by: 

E = imV 2 CJ] (1-1) 

If the density of the flowing air is p , then the kinetic energy per vo- 
lume of air, that has a velocity V , is: 

E = |pV 2 [Jm" 3 ] (1-2) 

If we consider an area A perpendicular to the wind direction (see fig. 1.1), 
then it may be seen that per second a volume A flows through this area. 
V is the undisturbed wind velocity. 




fig. 1.1 Area A and volume flow per second A. 

So the flow per second through A is: 

$ = V A [mV 1 ] (1-3) 

V oo 

The power that flows with the air, through area A, is the kinetic energy 
of the air that flows per second through A. 
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Pcwer = Energy per second. 



Power = Energy per volume * Volume per second. 



Equations (1-2) and (1-3) combined give: 



air 2 °° 



[Js 



-1 



= watt] 



* V A 



00 



P air* A 



[W] 



(1-4) 



This is the power available in the wind; as will be seen, only a part of this 
power can actually be extracted by a windmill. 

The above derived relation for the power in the wind (1-4) shows clearly that: 

- the power is proportional to the density p. 

This factor can not be influenced and varies slightly with the height 
and temperature (For 15°C at sea level p = 1,225 kgm 

- in case of horizontal axis windmills the power is proportional to the 

2 2 
area A = ttR (area swept by the blades) and thus to R . Radius R is cho- 
sen in the design. 

- the power varies with the cube of the undisturbed windvelocity V . Note 
that the power increases eightfold if the windspeed doubles. 

2. HORIZONTAL AXIS WINDMILL ROTOR-. 

To extract the power from the wind, several devices have been used and are 
still in use throughout the world. Examples of such devices are sailboats and 
windmills. 

This book deals only with the design of rotors for horizontal axis windmills, 
which are rotors with the axis of rotation in line with the wind velocity. The 
rotor rotates because forces are acting on the blades. These forces are acting 
on the blades because the blade changes the air velocity. The next paragraph 
deals with the relations between the velocity at the rotor blade and the for- 
ces acting on the blade. 



2.0 Airfoils. 



The rotor of a windmill consists of one or more blades attached to a hub. 
The cross sections of these blades can have several forms, as illustrated 
in fig. 2.1 and we call these cross sections airfoils. 







airflow ^ ^^^^^^^^^^ 

flat plate 


airflow ^_ ■ -<^^^^^^ 

arched plate 


symmetric airfoil 


airflow ^ ^ , '^ Vss s 1 
sail with pole 


m' rf 1 nw ^, ^^^^^^^^ 


airflow ^ ^ — 


cambered airfoil 


sail with pole 


airflow "^^^^Sv 


airflow fc ^l!^^^ 


highly cambered airfoil 

- 


sailwing 



fig. 2.1 Types of airfoils. 



An airfoil is a surface over which air flows. This flow results in two 
forces: LIFT and DRAG. Lift is the force measured perpendicular to the 
airflow - not to the airfoil I Drag is measured parallel to the flow. 
See fig. 2.2. h 



airflow 




fig. 2.2 Lift and drag. 



All airfoils require some angle with the airflow in order to produce 
lift. The more lift required, the larger the angle. 

The chord line (fig. 2.3) connects the leading edge and the trailing 
edge of the airfoil. The angle required for lift is called angle of attack ex. 
The angle of attack is measured between the chord line and the direction 
of the airflow. See fig. 2.3. 




fig. 2.3. Chord line and angle of attack. 



We want to describe the performance of an airfoil independent of size and 

2 

airflow velocity. Therefore we divide lift L and drag D by £pV A where 



p = air density 
V - flow velocity 

A = blade area (= chord * blade length) 



Ckgm 3 ] 
Cms ! ] 
Cm 2 ] 



The results of these divisions we call lift coefficient and drag coeffi- 
cient C, "* 



c l= 2 L 


> 


> 

[-] 


c - D 




[-] 


d *PV 2 A 







(2-1) 



(2-2) 



As stated before, the amount of lift and drag that is produced, depends on 
the angle of attack. This dependence is a given characteristic of an airfoil 
and is always presented in C^-a and c ^~ c ^ graphs. See fig 2.4. 



For some airfoils the chord line is defined otherwise. 
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fig 2.4. Lift and drag characteristics. 



In Appendix II C^-ci and C^-C^ characteristics of a series of NACA airfoils 

is presented as an example. For the design of a windmill it is important to 

find from such graphs the and a values that correspond with a minimum 

C./C -ratio^pThis is done in the following way: 
a 1 

In the C^/C^ graph a tangent is drawn through C^=C^=0. See fig 2.4.b. 

From the point where the tangent touches the curve, we find and . Fronr 

fig. 2. 4. a we find the corresponding angle of attack a. 

The and a values that are found in this way we call C^-design and a-design 

and the division of C. by C, is the minimum C ,/C, -ratio: (C./C,) . . 

d J 1 d 1 d 1 mm 

Table 2.1 on p. !0 gives these design values for several airfoils. 
Note that it is not important for the behaviour of the airfoil whether it is 
standing still in an airflow with velocity W or that it is moving with velo- 
city W in air that is at rest; what matters is the relative velocity that is 
"seen" by the airfoil. See fig. 2.5. 





fig 2.5 T-lative velocity. 



A blade element of a windmill rotor "sees" a relative velocity that results 
from the wind velocity in combination with the velocity wit-t vhich '„Pe hlade 
element moves itself. See fig 2,6. 



a;) We will see on p. 11 that a maximum power is obtained when the drag to lift 
ratio is as small as possible. 



air velocity 



Fig. 2.6 Relative velocity on rotor blade 



$ is the angle between the relative velocity W and the rotor plane, 



TABLE 2 . 1 



airfoil 
name 



geometrical 
description 



sail and pole 



c/10 



73 



^ -« — — 



flat steel plate 



<V c iii 



0.1 



0. 1 



0.8 



0.4 



arched steel plate 




f/c=0.07 
f/c=0. 1 



0.02 
0.02 



arched steel plate with 
tube on concave side 




d f/c=0.07 
d<0 ' lc f/c-0.1 



0.05 
0.05 



arched steel plate with 
tube on convex side 




f/c=0. 1 



0.2 



cloth or sail 



sail wing 




0.05 



tube steel cable 



sail trouser 

f/c=0.1 d 



tube 



=0.6f 



,c/4 cloth or 
sail 



0. 1 



0.9 
I .251 



5 0.9 
4 1.1 



14 



1 .25 



1.0 



.0 



NACA 4412 



see appendix II 



0.01 4 



0.8 



NACA 23015 



see Lit(I) in appendix I 



0.01 



0.8 



II 



1 Torque and pow e r characteristics 



L = Lift 



The components in the plane of rotation of the lift forces result in a 
force working in tangential direction at some distance from the rotor 
center. This force is diminished by the component of the drag in tan- 
gential direction. The result of these two components is a propelling 
force in tangential direction at some distance from the rotor center. 

The product of this propelling 
force and its corresponding 
distance to the rotor center 
is the contribution of the 
blade element under considera- 
tion to the torque Q of the 
rotor. The rotor rotates at an 
angular speed ft (= 2tt* number 
of revolutions per second). 



n = 2irn [rad sec l ] (2-3) 

The power that such a rotor ex- 
tracts from the wind is trans- 
formed into mechanical power. 
This power is equal to the 
product of the torque and the 
angular speed. 




component of 
drag in tangential 
direction 



'component of lift in 
tangential direction 

"relative velocity W 



Q - torque 
; fl = angular speed 



power P = Q * n 



[Nm] 

[rad s } 
CW] (2-4) 



A windmill of given dimensions transforms kinetic energy from the wind in- 
to a certain amount of power. Equation (2-4) clearly shows that a windmill 
for a high torque load (for example a piston pump) will have a low angular 
speed; a high speed design will only produce a small amount of torque (for 
example for a centrifugal pump or an electricity generator). 

We call a graph, that shows the dependance of the windmill torque on the angu 
lar speed, a windmill torque characteristic. Fig. 2. 7. a shows torque character 
tics of two different windmills designed for the same power but for different 
angular speeds. The windmill torque characteristic depends on wind speed V^, 
so we have many curves in one characteristic. 
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rotor radi : -Urn) 



Mima ) 



Q (Km) 



2 4 6 6 10 »il 

fig. 2. 7. a. low speed windmill 

torque characteristic 

p .(w) 

R= 1 (m) 




(rads ) 

high speed windmill 
torque characteristic 




V <ms~ l ) 




fig. 2.7.b. low speed windmill 
power characteristic 



high speed windmill 
power characteristic 



With relation (2-4) it is very simple to derive from the torque characteristics 
the corresponding power characteristics. See fig. 2.7.b, where the power-angular 
speed curves, belonging to windmills of fig. 2. 7. a, are shown. 
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Note: I) The power of the two windmills is the same, but is 
delivered at different angular speeds ft. 

2) The maximum power is delivered at a higher angular 
speed than the maximum torque. 

3) The maximums of the power curves in fig. 2.?.b. va- 
ry with the cube of the angular speed ft: 

P ~ft 3 (2-5a) 
max 

while the corresponding torque values vary with the 
square of the angular speed ft: 

Q (at P = P ) ~ " 2 (2 - 5b > 
max 

4) Tlie starting torque, i.e. the torque at zero revolutions 
per second, is considerably lower for high speed than 
for low speed windmills. 

Before selecting the speed of the rotor to be designed, the designer must 
compare the torque characteristic of the load with the torque characteristic 
of the rotor. 

For a. proper match of a load to a windmill rotor it is important that both 
load and windmill operate at angular speeds where their efficiencies are maxi- 
mum. The angular speed at which the rotor has its maximum efficiency is not 
always equal to the angular speed at which the load has its maximum efficien- 
cy. In that case we need a transmission. It is in most cases not difficult 
to determine the transmission factor needed for the optimum angular speeds 
of both load and rotor, but remember that a transmission also changes the tor- 
que. In practice this means that the transmission factor cannot be chosen on 
the basis of angular speeds only. 
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2,2 Dimens ionless coefficients. 

In order to be aDle to compare the properties and characteristics of 
different windmill designs under different wind conditions we write 
the mechanical power as the power in the air multiplied by a factor C 



P . = C * P . 
mech p air 



(2-6a) 



0^ is called power coefficient and is a measure for the success we have 
in extracting power from the wind. With relation (1-4) we may write 

(2-6b) 




For the same reasons we divide the speed u of the rotor at radius r by 
the windspeed. See fig. 2.8 

L R 




u= wQr 



fig. 2.8 Definition of speed ratio. 



The 



result Oy ) we call local speed ratio and is noted 



u Sir 
r " V ~ V~ 

CO CO 



(2-7) 



The speed ratio of the element of the rotor blade at radius R we call 
tip- speed ratio: 



(2-8) 



Note: Later it will be shown that a windmill has one valup of > at which 
the power coefficient is maximum. This X is often called 'the tip-speed 
ratio of a windmill 1 or 'the speed ratio of a windmill'. 
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There is of course a direct relation between X and X^. Relations (2-7) 

and (2-8) together give 

A = I * X (2-9) 

r R 



From relation (2-4) we know *:hat 



« - 



(2-10) 



With this relation we define a d,imeris ion less torque coefficient in the 
following way: 



P = C JpV 3 ttR 2 



xv 

a 



Q = n 



(2-6b) 

(2-8) ^ 
(2-10) 



X JpV^ ttR^ 



We define: 




(2-11) 



Note that in this way relation (2-4) is still valid but now in dimen- 
sionless form: 



C = C * A 
P Q 



(2-12) 
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2.3 Basic form of a windmill characteristic. 

The power coefficient C p in equation (2-6) is not an efficiency but 
may be interpreted as a measure of the success that a windmill has in 
transforming wind energy into mechanical energy. For one specific wind 
mill Gp varies with the tip-speed ratio of the windmill. In dimension- 
less form this is shown in a so-called Cp - X characteristic based on 
formulas (2-6) and (2-8). 

See fig. 2.9 where one curve now represents all the curves for dif- 
ferent V of fig. 2.7.b. 




fig 2.9 C - X characteristic. 



This characteristic is independent of air density p, windspeed V m and 
radius R. 

Using relation (2-12) we may derive from fig 2.9 a dimensionless form 
of the torque-speed characteristic of the windmill: C^-X curve; see 
fig 2.10. Also here one curve represents all curves of fig 2. 7. a. 




I 2 J i 5 6 7— A 



fig 2.10 Cq-X characteristic. 

Note that the power is zero if X * but that the torque is not. 
See relations (2-4) and (2- J 2). 



17 



2.4. Maximum power coefficient. 

The power coefficient Cp as defined with relation (2-6) describes 
how much power we get from the wind with a windmill. The power in 
the wind is given by relation (1-4). We are of course very interes- 
ted in how much wind power we can transform into mechanical power 
with a windmill. In other words, we want to know what the highest 
power coefficient is for a given windmill that is designed for 
a certain tip-speed ratio. Betz was the first one to show that the 
theoretically maximum attainable power coefficient is 0.593. This 
result will be clarified in the next paragraph. Three other effects 
cau*e a further reduction of the maximum power coefficient. How to 
find the maximum power coefficient that takes these effects into 
account will be explained in the next four paragraphs: 



2^4.0 Betz coefficient 

2 J.. I Effect of wake rotation on maximum power coefficient 

2 .".2 Effect of C,/C, -ratio on maximum power coefficient 
d L 

2.4.3 Effect of number of blades on maximum power coefficient 



1 
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2.4.0 Betz coefficient 

It is not possible to transform all the wind energy that flows 
throughcross sectional area A (fig 1.1) into mechanical energy. 
If we could transform all the energy in the air this would mean 
that we could extract all kinetic energy from the air; the air ve- 
locity behind the rotor would then be zero and no more air would 
flow through the rotor. The process of extracting kinetic energy 
from the wind will stop and no more power will be transformed. If 
on the other hand the air velocity behind the rotor is equal to the 
wind velocity, no kinetic energy has been extracted and also in this 
case no power will be transformed. In this way it may be understood 
that, if the flow velocity behind the rotor is either zero or equal 
to the wind velocity V , in both cases the mechanical power is zero. 
Between these values there is an optimum value of the wind velocity 
behind the rotor. 

Betz found this value to be ^ and calculated the maximum power 
coefficient (Betz coefficient). 

C p =^=0.593 (2-13). 
max 

This value is however only valid for a theoretical design for a high 
tip-speed ratio , with an infinite number of blades and a blade drag 
equal to zero. The effect of deviations from these three assumptions 
will be shown one by one in the next three paragraphs. 
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2.4.1 Effect of wake rotation on maximum power coefficient. 

The Betz coefficient suggests that, independent of the design tip-speed 
ratio X q , we may expect a maximum power coefficient C p of 0.593. 
Relation (2-13) is however only valid for high tip-speed ratios and 
for low tip-speed ratios considerable deviations exist. This can be 
explained in the following way: 

The power is: torque * angular speed. The torque is produced by for- 
ces acting on the blades in tangential direction, multiplied by their 
corresponding distances to the rotor center. These forces are the re- 
sult of velocity changes of the air in tangential direction (action = 
reaction; force = mass * velocity change pei: unit of time). The direc- 
tion of the velocity change in the air is opposed to the direction of 
the forces acting on the blades. Since the air has no tangential velo- 
city before passing the rotor, the velocity change means that behind 
the rotor the wake rotates in a direction opposite to that of the rotor. 
This wake rotation means a loss of energy because the rotating air con- 
tains kinetic energy (see relation 1-1). Since a certain amount of po- 
wer is to be transformed, we know from relation (2-4) that a low tip- 
speed ratio (= low angular speed ft), means that the torque Q must be 
high. High torque means large tangential velocities in the wake; the 
consequence is a loss of energy and a lower power; the more so if the 
design tip-speed ratio is lower. The result is shown in fig, 2.11. 
The graph presented in fig. 2.11 shows the collection of maximum ob- 
tainable power coefficients of ideal windmills i.e. wir 'lis with an 
infinite number of blades without drag. 




fig. 2.11 Collection of maximum power coefficients of ideal 
windmills . 
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2.4.2 Effect of C./C -ratio on maximum power coefficient, 
a 1 

The factor C,, as defined in par. 2.0, is a measure for the resistance 
d 

of the blades against moving through the air. The C^/C^-ratio deter- 
mines the losses due to this resistance. These losses are calculated 
and included in the collection of maximum power coefficients of fig. 
2.11. The results are shown in fig. 2.12. 




X 



'. 2 3 4 5 6 7 8 9 10 1 1 12 13 ! 4 15 



fig. 2.12 Effect of C /Cj-ratio on C^-max for 

a rotor with an infinite number of blades. 



Fig. 2.12 shows that a rotor, designed for a tip-speed ratio A = 2.5 
with airfoils having, for example, a minimum C^/C^-ratio of 0.05, will 
have a maximum power coefficient Cp = 0.46. If the rotor is designed 
for A = 10, at the same C^/C^ value, the value at A = 10 will have 
a maximum of 0.3. 

Note that from fig. 2.12 it is clear that it is useless to design a 

rotor for X = 10 with airfoils that have (C./C.) . =0.1. 

d 1 mm 
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2.4.3 Effect of number of blades on maximum power coefficient. 



The numliar of blades also affects the maximum power coefficient. 

This is caused by the so-called "tip-losses" that occur at the tips of 

the blades. These losses depend on the number of blades and the tip-speed 

ratio. The losses have been calculated and as example are included in the 

collection of maximum power coefficients for (C./Cj . = 0.03 of fig. 

d 1 mm 

2.12. The results are shown in fig. 2.13. 




fig. 2.13 Influence of number of blades B on C -max for C,/C = 0,03 

r d 1 

In appendix III graphs like fig. 2.13 are shown giving the maximum expec- 
ted power coefficient for A-design between 1 and 15, 

The first group of graphs shows C for constant number of blades B 

max 

while the C,/C -ratio is varied. The second group gives C for constant 

, . max 

C./C -ratio while B is varied, 
d 1 

Conclusion: If design figures for tip-speed ratio X, number of blades B 
and C^/C^-ratio have been chosen, the expected power coefficient C,, may be 
read from the graphs in appendix III. 
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3 DESIGN OF A WINDMILL ROTOR. 

3.1' Calculation of blade chords and blade setting. 

In chapter 2 it was shown that the selection of the number of blades B af- 
fects the power coefficient. Although B has no influence on the tip-speed 
ratio of a certain windmill, for the lower design tip-speed ratios, in ge- 
neral, a higher number of blades is chosen (see table 3.1). This is done 
because the influence of B on C is larger at lower tip-speed ratios. A 

P 

second reason is that choice of a high number of blades E for a high design 
tip-speed ratio will lead to very small and thin blades which results in 
manufacturing problems and a negative influence on the lift and drag properties 
of the blades (this problem will be dealt with in chapter 4). 



X 

o 


B 


1 


6-20 


2 


4-12 


3 


3-6 


4 


2-4 


5-8 


2-3 


8-15 


1-2 



table 3.1 Selection of number of blades. 

A second important factor that affects the power coefficient is the drag. 

Drag affects the expected power coefficient via the C ,/C, -ratio . This will 

d 1 

influence the size and, even more, the speed ratio of the design. 
In paragraph 2.0 a list is shown with several airfoil types md their cor- 
responding minimum C^/C^-ratios (table 2.1 on p. 10). Promising airfoils in 
this table have a minimum C/C -ratio between 0.1 and 0.01. 
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A large C^/C^-ratio restricts the design tip-speed ratio. At lower tip-speed 

ratios the use of more blades compensates the power loss due to drag. See 

Appendix III. In this collection of maximum power coefficients it is seen 

that for a range of design speeds 1 <A Q < 1 the maximum theoretically attainable 

power coefficients lie between 0.35<C_, ^0.5. 
r *max 

Due to deviations, however, of the ideal geometry and hub losses for example, 
these maximums will lie between 0.3 and 0.4. This result shows that the choice 
of the design tip-speed ratio hardly effects the power output. Two other fac- 
tors, however, limit the choice of the design tip-speed ratio. One is the cha- 
racter of the load. If it is a piston pump, scoop wheel or some other slow 
running load, that in most cases will require a hi^h starting torque, the de- 
sign speed of the rotor will usually be chosen low; this allows the designer 
to use simple airfoils like sails or steelplates. If the load is running fast 
like a generator or a centrifugal pump, then a high 'sign speed will be se- 
lected and airfoils with a low C^/C^-ratio will be preferred. The second fac- 
tor is that the locally available technologies will often restrict the possi- 
bilities of manufacturing blades wi.f.h airfoils having low C^/C^-ratios . But 
even in the case of a high speed design, simple airfoils like arched steel pla- 
tes, can give very good results. 

We will pass over in silence the problems related to starting 
torque and optimum angular frequency of the load. Now we can. design a windmill 
rotor for a given windspeed V and a power demand P. Selecting in table 2.1 an 
airfoil in terms of a minimum expected C^/C^-ratio, we can choose a design 
tip-speed ratio \ Q with the help of appendix III. With table 3.1 a number of 
blades B is chosen and, returning to appendix III we can find the maximum power 

coefficient C that may be expected, 

^max 

example: design rotor with 7% arched steelplates. 

Table 2.1 - (0^)^ = 0.62 - 

Appendix III Page 48 - 1<X Q <#; we choose \ q = 4. 

Table 3.1 - B - 4 

Appendix III Page 43 - C p =0.48 

max 
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Willi relation (2-6) we can now calculate Uic desired radius R of the rotor. 
For conservative design we take C = 0.8 * C 

max " 









/ " 











example: the rotor to be designed mi.st deliver 2100 Watts of mecha- 
ni-jal power CP - HQQViutt) in a wind with V - 8 ims ]. 



/? - / — 1.7 \ m J 



Design of the blades. 

We need the following data: 



rotor radius R [ml 

number of blades B [-■] 

design tip-speed ratio \ [~ ] 
airfoil data: 

design lift coefficient C [- 1 

o 

corresponding angle 

of attack a_ [-] 



Airfoil data may be found in table 2.1, Appendix II and literature (1), 
(5), (f>) ar.d (12). 

Onr - Lhese data are known, it is now very simple to calculate the 
blade geometry; i.e. the chord c of the blade and the blade angle 3, 
the angle between the chord and the plane of rotation, fig. 3.1. 




W 

fig. 3.1 Made setting 3. 



Only three simple formulas are needed and one graph 

X r = X o X r/R (2 ' 9 > 

8TTr ,i ^ 
C = BC~ ( " COS * J < 3 * 

6 = * - a (3.2) 
and graph - <J» (appendix V) 

The underlying theory is too complicated to be explained here (see 
appendix IV and the literature references appendix I), The reader who 
is primarily interested in the design of the rotor, can do without 
this theorv. 



Kow the design procedure is as follows: 

Divide the blade with radius R in a number of parts of equal length. 

In this way we find cross sections of the blade. Each cross section 

has a distance r to the rotor center and has a local speed ratio A » 

r 

according to (2.9). 

In appendix V we can find the corresponding angle $ (see fig. 3.1) for 
each cross-section. 4» is the angle of the relative air velocity W that 
meets the blade section at radius r. We now calculate the chord with 
relation (3.1). (For ease, (1 - cos <fc) has been added in the graph of 
appendix V). The blade angle at the corresponding radius is found 
with (3.2). 

example: We continue our design of a rotor with 

R - 1.7 Irnl 

B - 4 

X = 4 
o 

airfoil = 7% arched steel plate 
table 2.1 gives for this airfoil: 

= 0.9 (i.e. the value for minimum 

o 

a Q - 4° (i.e. the corresponding angle of attack). 

With equations (2-9) 3 and (3-1), (3-2) we can now com- 
pute the values of the following table: 
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Table 3.2 



cross 
section 
number 


r(m) 


X 

r 


4° 




ot 

o 


e° 




I 


0.2125 


0.5 


42.3 


4 


38.3 


0.386 


2 


0.4250 


1 


30.0 


4 


26.0 


0.398 


3 


0.6375 


1.5 


22.5 


4 


18.5 


0.399 


4 


0.8500 


2 


17.7 


4 


13.7 


0.281 


5 


1 .0625 


2.5 


14.5 


4 


10.5 


0.236 


6 


1.2750 


3 


12.3 


4 


8.3 


0.204 


7 


1.4875 


3.5 


10.6 


4 


6.6 


0. 177 


8 


I .7000 


4 


9,4 


4 


5,4 


0. I 59 



The result is the blade chord c and *\e blade setting g at various stations 
along the blades. Plotting the chords (fig 3.2) gives the blade form and 
plotting g shows the desired twist of the blade. 



rotor center 



cross section 
number 



/ l ; 

v c 

r 


2 3 4 5_ _6 7 

'_ - " 

R 







30 l 



20 



10' 



— I 

0.5R 



r (m) 
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3.1 Deviations from the calculated chorda and blade sett ing. 

In the last paragraph we showed how to calculate the ideal blade form. 
The chords as well as the blade angles as calculated in par. 3.0 vary 
in a non-linear manner along the blade. Such blades are usually diffi- 
cult to manufacture and lead to an uneconomic use of materials. In order 
to reduce these problems it is possible to linearize the chords and the 
blade angles. This results in a small loss of power. If the linearization 
is done in a sensible way the loss is only a few percent. 
In considering such linearizations it must be realized that about 75% 
of the power that is extracted by the rotor from the wind, is extracted 
by the outer half of the blades. This is because the blade swept area 
varies with the square of the radius; also the efficiency of the blades 
is less at small radii, where the speed ratio is small. On the other 
hand, at the tip of the blade the efficiency is low, due the so-called 
tip losses discussed in par. 2.4.3. 

For the reasons mentioned above, it is advised to linearize the chords c 
and the blade angles 3 between r = 0.5R and r = 0.9R. 



example: 



we linearize the blade chords a and angles as cal- 
culated in table 3.2 of par. 3.0. 



The nearest value of r to 0.5 R in table 3.2 is 
r = 0.85 (= 0,5 R). 

The nearest value of r to 0.9 R in table 3.2 is 
r = 1.4875 (= 0.875 R) . 

In table 3.2 for these values of r the following va- 
lues of c and were calculated: 



0.85 



0.281 13.7 
0.159 6.6° 



1.4875 
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We can now linearize the chords and blade angles by writing 
c and 8 in the following way: 

c = a 2 r + a 2 

S - a 2 r + a 4 

With the values of c and 6 at r = 0.85 xnd r = 1.4875* the 
constants a^ 3 a^ 3 a^ and a^ are found. 

c = -0.191 r + 0.444 
S - -11.14 r + 23.17 

Suppose we have a hub for the rotor and the hub has a radius 
r = 0.17 Cm] then we can calculate the chords and blade angles 
at the foot and at the tip of the blade: 

c„ + = -0.191 * 0.17 + 0.444 = 0.412 Cm] 
foot 

c,. - -0.191 * 1.7 + 0.444 = 0.119 Cm] 
tvp 

8- . = -11.14 * 0.17 + 23.17 = 21. 2° 
foot 

ft . = -11.14 * 1.7 + 23.17 - 4.3° 
tvp 



The result of the linearizations is shown in fig. 3.3 where blade form and 
twist in the blade are shown and compared with the blade form of fig. 3.2, 




fig. 3.3. Linearized chords and twist. • r(m) 
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As may be seen, from fig. 3.3 the changes in chords and blade angles are 
very small at the outer half of the blades. At the inner half the chords 
also remain almost unchanged. A rather large change is found in the blade 
angles for r < 0.5 R. For reasons as stated on the first page of this chap 
ter, this will not lead to any significant power loss but may have a consi 
derable effect on the torque that is produced at low angular speeds , In 
general the starting torque will be less and in cases where the starting 
capacities of the windmill are very important, this effect should not he 
forgotten. An example of a load that demands a high starting torque is 
the single stroke piston pump. For this load, the size of the windmill is 
often determined by the demanded starting torque. 
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4. EFFECT OF THE REYNOLDS-NUMBER. 

4 .0 Uependance of airfoil characteristics on the Re-number. 

The airfoil characteristics depend on the so-called Reynolds'-number (Re) 

of the flow around the airfoil. For an airfoil Re is defined as Re = Hl£ , 

v 

where W is the relative velocity to the airfoil, c is the chord and v is 
the kinematic viscosity (in our case that of air). 

All airfoils have a critical Re-number. If the Re-number of the flow a- 
round the airfoil is less then this Re cr £ t i cal then the C-^-value is lower 

and the C, -value is higher; above this Re . . the performance is con- 
u critical 

siderably better. See for example fig. 4.1 where the effect of the Re-num- 
ber on (C ,/C. ) . is shown, 
d 1 mm 




\ — i — i — i — , — i — i — i — i — i 

tO 90 SO 80 tOO T20 110 160 TSO H^Re 



fig. 4.1. Effect of Re-number on (C./C.) . -ratio for three dif- 

d l mm 

ferent airfoils. 

In general the critical Re-number for airfoils with a sharp nose will 
4 

be 10 while for the more conventional airfoils like NACA the critical 
Re-number is about 10^; some of the very modern airfoil types have a 
critical Re-number of about 10^. 

Fig. 4.2 shows the inverse value of the C./C -ratio of various airfoils 

d 1 

as f (Re). 
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fig. 4.2. Inverse value of minimum C J /C -ratio as function 

d. 1 

of the Re-number for several airfoils. 
From lit(5). 
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4.i Calculation of the Re-number for the blades of a windmill rotor. 



For the condition that the rotor runs at A=X . the Re-number of 

optimum 

the flow around the airfoil can be determined with fig. 4. 3 in the fol- 
lowing way: 

if B = number of blades 

r ~ radius = distance to rotor center of blade element under 
consideration 

A r = speed ratio of blade element under consideration 

= design lift coefficient of blade element under consi- 
deration 
= undisturbed windspeed 

the Re-number is: 

V a * r 

Re = — * Re 

B * C N 

Re^ may be read from the graph presented in fig. 4.3 (valid 

— f> 7 — \ 

for air: kinematic viscosity „■ = 15 * 10 [m s ]) 
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example: we will check the Re-number for the rotor as designed in 
chapter 3. 

airfoil: 7% arched stesl plate 
radius : R - 1.7 [m] 
tip-speed patio: \^ - 4 
C j- design = 0.9 
Number of blades B = 4. 



1) at the tip X - X - 4. 

r = R = 1.7 M 
4 



fig 4. 3- R& R - 9*10 



Re r=R = Re N * 



R * 7 9 * 10 4 * 1.7 



B * C l 4 * 0.3 



Re D - 4.25 *10 4 V 
r=R ° 



2) at r = 0.5 R = 2; r = 0.85 



17 * 10 4 * 0.85 
Re _ n CD - V = 4 * 10 V 

r=0 ,5R *nn °° ° 

4 * 0.9 



3) at r = 0.2 R \, T = 0.8 r ~ 0.34 



28 * 10 4 * 0.34 

'r=0.2R ~ 



Re n nn - - 2.6 * 10 4 V 



4 * 0.9 

Conclusion: for the whole blade 3 the Re^numler is 3 even 
for very low windspeeds 3 higher than the critical Re-num- 
ber for steel plates (= 10 4 ) . 

Thus the assumed minimum C^/C^-ratio is correct . 
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KACA 4412 



(Stations j.jiu uiniiwti's gut'ii m 
per cent of airfoil chord) 



Upper surface 


Luu cr 

- - 


surface 

_- .. 


Station 


UrdiiLiLf 


Elation 


Ordinate 














1.25 


2.44 


1.25 


- 1.43 


2.5 


3.39 


2.5 


- 1.95 


5.0 


4.73 


5.0 


- 2.49 


7.5 


r 

O.iO 


7.5 


— 2.74 


10 


6.59 


10 


- 2.S0 


15 


7-S9 


15 


- 2.SS 


20 


8.S0 


20 


- 2.74 


25 


9.41 


25 


- 2.50 


30 


9.76 


30 


- 2.26 


40 


9.S0 


10 


ISO 


50 


9.19 


50 


- 1.40 


i'tC 


S/.4 


GO 


- 1.00 


70 


e.69 


70 


- 0.65 


Si) 


4.S9 


SO 


- 0.39 


90 


2.71 


90 


- 0.22 


95 


1.47 


95 


- 0.1G 


100 


(0.13) 


100 


(- 0.13) 


100 





100 





L.E. radius: 1.58 






Slope of radius through L.E.: 0.20 



NAC.V 4415 



(Stations and orduwtes p;v>i 'm 
per rent (if :i ii full ( h.ii d) 



Upper 

.. 


surfarc 


Lower Kurface 


Station 


Odimite 


Station 


Ordinate 















1.25 


3 07 


1.25 


- 1.79 


2.5 


4.17 


2.5 


- 2.48 


5.0 


5.74 


5.0 


- 3.27 


7.5 


0.91 




— 6.1 1 


10 


7.8-1 


10 


- 3.98 


15 


9.27 


15 


■■ '.I; 


20 


10.25 




— 4.15 


25 


10.92 


25 


- 3.9S 


30 


11.25 


30 


— 3.75 


in 


11.25 


40 


- 3.25 


50 


10.53 


50 


•) j 2 


GO 


9.30 


GO 


- 2.14 


70 


7.G3 


70 


- 1.55 


80 


5.55 


SO 


- 1.03 


90 


3.08 


90 


- 0.57 


95 


1.67 


95 


- 0.36 


100 


(0.16) 


100 


(- 0.16) 


100 




100 





L.E. radius: 2.48 






Slope of radius through L.Li.: 0.20 



NAOA 4418 



(Stations and ordinate.* given in 
jut cent of nirfoil chord I 



Upper surface 


Lower surface 


Station 


Ordinate 


Station 


Ordinate 













1.25 


3.76 


1.25 


- 2.11 


2.5 


5.00 


2.5 


- 2.99 


5.0 


6.75 


5.0 


- 4.06 


7.5 


6.U6 


7.5 


- 4.67 


10 


9.11 


10 


- 5.06 


15 


J.U.00 


15 


- 5.19 


20 


11.72 


20 


- 5.56 


25 


12.10 


25 


O.-t J 


30 


12.76 


30 


- 5.26 


10 


12.70 


40 


- i.70 


50 


11.85 


50 


- 4.02 


60 


10.44 


60 


- 3.24 


70 


8.55 


70 


- 2.45 


80 


6.22 


80 


- 1.67 


90 


3.46 


90 


- 0.93 


95 


1.S9 


95 


- 0.55 


100 


(0.19) 


100 


f- 0.19) 


100 




1UU 





L.E. radius: 3.50 






Slope of ! udiu> tl 
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XACA 4421 
(Stations and ordinates given in 
per cent of airfoil chord) 



Upper surface 


Lower 


surface 


Station 


Ordinate 


Station 


Ordinate 













L25 


4 »5 


1.25 


— 2.42 


2.5 


Q.D-k 


2.5 


- 3.48 


5.0 


7.82 


5.0 


- 4.78 


7.5 


9.24 


7.5 


- 5.62 


10 


!0.35 


10 


- 6.15 


15 


12.04 


15 


— 6.75 


20 


13.17 


20 


- 6.98 




13.88 


25 


- o.W 


30 


14.27 


30 


- 6.76 


40 


14.16 


10 


- 6.16 


50 


13.18 


50 


- 5.34 


60 


11.60 


60 


- 4.40 


70 


9.50 


70 


- 3.35 


80 


6.91 


80 


- 2.31 


90 


3.85 




- 1.27 


95 


2 11 


r t.- 


- 0.71 


100 


f0.22) 


100 


(- 0.22; 


100 




100 





L.E. radius: 4.85 






Slope of radius through L.1 


:.: 0.20 



NACA 4424 
(Station- rmd ordinates given in 
per cent of airfoil chord) 



Upper surface 


Lower surface- 


Station 


Ordinate 


Station 


Ordinate 














0.530 


3.964 


1.970 


- u. 472 


1.536 


5.624 


3.464 


- 4.656 


3.775 


7.942 


6.225 


- 6.066 


6.153 


9.651 


8.847 


- 0,931 


S.Gil 


11.012 


11.389 


- 7.5 1 2 


13.674 


13.045 


16.326 


- S.169 


18.858 


14.416 


21.142 


- 8.41 G. 


21.111 


15.2.S7 


25.SS9 


- S. Ill 


29.401 


15.73S 


30.599 


- S.238 


10.000 


15.606 


10.000 


- 7.006 


50.235 


14.474 


49.705 


- 6.698 


60.405 


12.671 


59.595 


- 5.562 


70.487 


10.312 


09.513 


- 4.312 


80.164 


7.447 


79.536 


- 3.003 


90.320 


4.099 


S9.GS0 


— 1.1 to.) 


95. 190 


2. J !') 


91. Sill 


- OHM 


100.000 





11/0.000 


'.1 


L.E. radius: 6.33 






S!oj>r> of radius tl 


ironph L.' 


•:.: 0.20 
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APPENDIX IV 



Note on the theoretical assumptions on which the design method 
is based. 

The design procedure presented in this publication is based on general 
momentum theory and blade element theory as can be found in for example 
lit (2-3-6-10). This simple design procedure ignores tip effects as 
mentioned in paragraph 2.4.3. In the collection of maximum power coef- 
ficients given in AppendixIII, however, the tip losses are included. 



As described in paragraph 2.4.0 - 2.4.3, the attainable power coefficient 
can be described with the following effects: 

1) Betz coefficient C p - — 

Betz 27 

2) wake rotation 

3) blade drag 

4) tiplosses due to finite number of blades 

Effects 1) and 2) can be described with the following approximation 
(error < 0.5 percent for X > 1): 

-1 .29 

r - 16 _ e "0.35X 

ideal 

with X=X 

optimum. 

Effect 3) can be described by reducing C with an approximation for 

*ideal 

power loss due to drag (max error = 2 percent for X=l; errors 0.1 percent 
for X > 2.5) 

16 C d 

C P , = 27 C7 X 
drag 1 

C d C d 
with X=X . . and yr- is 77- at 

optimum C l C t ^gign 



Cp including drag is: 
max 



-1 1Q C 

r - r - c - 'e-°' 35x - II a^> 

S " C P., , V " 27 v C. A) 

max idea drag 1 

Effect 4) can be included by multiplying C with factor n : 

1 max a 



( . 1.386 „. <K2 



B = number of blades; $ is found in Appendix V with \ ^ - * t £ p 



Thus C p is 
max 



c p = n B (c p - c p ) 
max ideal drag 
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C p = (1 - — g— sin 7 ) 27 (e C7 X) 

max 1 



APPENDIX V 



\ - si ntp(?co3(fr-l ) 
r ( i-COS (jj) (2coo ^+1 



\ 

r* 
1 




( 


-cos$ ) 


10.00 


'.807 





00220 


9.75 


3.904 





00232 


9. 50 


4 . 006 





00244 


'.25 


4.113 





00258 


9.00 


4.2?? 





00272 


B. 75 


4. M7 


0. 


00288 


8. 50 


4.475 





00305 


6.25 


4.607 





00323 . 


8.00 


4.750 





00343 


7.75 


4. 902 





00368 


7. 50 


5.06 3 





00390 


7.25 


5.236 





0041 7 


7 . 00 


5.420 





00447 


b.75 


5.618 





00480 


6. 50 


S.831 





00517 


o.25 


6. 060 





00559 


6.00 


fa. 308 





00605 


5.75 


6.577 





,00658 


5-50 


6.870 





00718 


^.25 


7. 190 





,00786 


5.00 


7.540 





00865 


4.75 


7.926 





,00955 


4.50 


8. 353 





,01061 


4.2S 


6.8?7 





,01 184 


4.00 


9.357 





01336 


5.75 


9.954 





,01505 


3-50 


10.630 





,01716 


3.25 


11.402 





,01774 


3. 00 


12.290 





,02292 


2-75 


13.322 





,02691 


2.50 


14.534 





,03200 


2.25 


15.975 





,03862 


2.00 


17.710 





,04739 


1.75 


19.830 





,05930 


1,50 


22.460 





,07585 


1.25 


25.773 





,09948 


1,00 


30.000 





,13397 


75 


35.420 





,18507 


0.50 


42.290 





,26025 


0.25 


50.642 





,36584 


0.00 


60.000 





,5 












Angle $ between relative velocity and plane of the rotor versus 
the speed ratio of an element at radius r for a windmill with a 
flow equal to the flow of an ideal windmill. 
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